o 
o 

(N 



c3 



Fractals, Vol. 0, No. (0000) 

© World Scientific Publishing Company 



FAST RELAXATION TIME IN A SPIN MODEL 
WITH GLASSY BEHAVIOR 

GIANCARLO FRANZESE 

Center for Polymer Studies and Department of Physics, Boston University 
590 Commonwealth Avenue, Boston, MA 02215, USA 



Abstract 

We consider a frustrated spin model with a glassy dynamics characterized by a 



a 

fj ■ slow component and a fast component in the relaxation process. The slow process 

involves variables with critical behavior at finite temperature T p and has a global 
character like the (structural) a-relaxation of glasses. The fast process has a more 
local character and can be associated to the /3-relaxation of glasses. At temperature 
T > T p the fast relaxation follows the non-Arrhenius behavior of the slow variables. 

ON ! At T < T p the fast variables have an Arrhenius behavior, resembling the a — (3 

bifurcation of fragile glasses. The model allows us to analyze the relation between 

L_. the dynamics and the thermodynamics. 



a 

O 
O 



1. INTRODUCTION 



A recent reviewEJ has pointed out some of the leading questions in the theoretical study of 
glassforming liquids. Following the fruitful path of theoretical studies on spin models with 
glassy behaviora, here we present results on a frustrated spin model that give an insight 
on some of these questions, showing relations between the dynamics and the statics of the 
model. 

We consider a model, introduced by Coniglio and coworkersauBHoH, with two kinds of 
variables, coupled to each other. One of them has a finite temperature transition and a 
slow relaxation, while the other has a transition only at zero temperature in two dimensions 
(2D) and has a fast relaxation at finite temperature. The model has been showM to have 
glassy behavior and the slow and fast components of the relaxation can be related to the 
a and (3 processes, respectively, of glasses. Here we show numerical results resembling the 
a-f3 bifurcation processa observed in experiments on fragile liquidsa'Ej. 

In the following we give an introduction to the problem. In Sec. 2 we introduce the spin 
model, reviewing some results on its phase diagram and on its glassy behavior. In Sec. 3 
we show some numerical results about the bifurcation in this model and the connection 
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with the thermodynamics. In Sec. 4 we give the conclusions and the perspectives for future 
work. 



1.1 Glasses And The a- (3 Bifurcation 

The majority of liquids can form a glass if cooled at high enough rate. The glass state is 
experimentally defined by measuring macroscopic quantities, like the diffusion coefficient or 
the viscosity. The calorimetric glass transition temperature T g is defined as the temperature 
at which the viscosity reaches a value of 10 12 sPa and depends on the cooling rate. 

The glass transition is considered as a dynamic off-equilibrium phenomenon, but a large 
amount of work has recently shown that the dynamics of glassy systems is strictly related 
to its static properties, i.e. to its thermodynamics!!!. The main focus of these studies is the 
slow component of the dynamics. One of the aim of this paper is to show that also the fast 
processes can be related to the thermodynamics of the system. 

Following the Angell's classification^, glassforming liquids are divided in 'strong' and 
'fragile'. Strong liquids are characterized by a relaxation time t(T) that increases with 
decreasing temperature T in a way well described by an Arrhenius law 

T(T)=T 00 exp[A/(k B T)] (1) 

where Too is interpreted as a characteristic microscopic relaxation time at infinite T and 
A is an activation energy (energy barrier) for global rearrangements (k B is the Boltzmann 
constant). A strong glassformer is, for example, SiC>2. 

Fragile liquids are defined as those whose t(T) shows a large departure from the Arrhenius 
law (non- Arrhenius behavior) for T > T g . For T <T g the relaxation time is well described 
by an Arrhenius law. Example of fragile glassformers are glycerol and orto-ther phenyl. A 
function widely used to fit data for fragile liquids is the Vogel-Tamman-Fulcher (VTF) law 

r(r) = r oo exp[B/(T-r )] (2) 

where r^ is a hight-T characteristic time, B and To are fitting constants. Since To < T g 
the diverging behavior is never reached. 

An alternative non- Arrhenius behavior is derived in the Mode Coupling Theory (MCT)I12I 
for T > T g , in which the relaxation time increases as a power law of the temperature, 

t(T) ~ (T - T MCT )- y , (3) 

where Tmct is the MCT transition temperature (with Tmct < T g ) and 7 is a parameter. 
The non- Arrhenius behavior can be explained, in a consistent way with the MCTtH, also 
by the Adam-Gibbs theory of the excluded volume, in which the relaxation time is 

r = r^ exp(C/TS c ) (4) 

where t^ ~ 10~ 14 s, C is a constant and S c is the entropy difference between liquid and 
crystal - 'excess' entropy. Since S c is, usually, constant below T g , the low-T Arrhenius 
behavior is recovered. 

This phenomenology concerns the structural a-relaxation, i.e. the process described by 
the long-time part of the relaxation functions of the system. For high T, the relaxation 
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functions, like the density-density correlation function, have a simple exponential behavior 
('simple liquid' regime) 

/(t) = / exp(t/r), (5) 

with /o and r depending on T and where t is the time. 

Decreasing T the relaxation functions usually show a 'two-steps' behavior where a fast 
relaxation process, associated to the first step, is followed by a slow relaxation process, the 
second step. The second step is the a-relaxation and is associated to a macroscopic process 
involving non-local degrees of freedom (as the global reorganization of the system p. Below 
a characteristic temperature T* > T g , the second step is well approximated in most cases 
by the Kohlrausch- Williams- Watts (KWW) stretched exponential function 

fKWw(t) = foe X p[{t/Tf] , (6) 

where /o, r and (3 depend on T. The KWW function recovers the simple exponential for 
(3 = 1, with the parameter (3 describing the departure from exponentiality. The fit of 
glassy relaxation functions, usually, shows (3 decreasing with the temperature, with (3 < 1 
for T < T*. In many cases, a better approximation of the second step of the relaxation 
function is given by using the OgielskiEJ stretched exponential function 

fo(t) = t- x f KWW (t), (7) 

with x depending on T. 

The first step (or /3-relaxation) is usually associated to any microscopic process^ occurring 
at very short time scale. In the MCT for atomic liquids it is associated to the fast diffusion 
of mobile particles inside 'cages' of immobile particleal^l. In molecular liquids it can be 
associated also to the relaxation of internal degrees of freedomEJ. 

The /3-relaxation presents some non-universal features in fragile liquidstl, and sometimes 
in intermediate liquidate, that are not well understood. The /3-relaxation time follows the 
non-Arrhenius a-relaxation behavior at high T and shows a crossover (the a — (3 bifurcation) 
to an Arrhenius behavior at a temperature T a -p ~ Tmct°- To have an insight on this 
phenomenon and its connection with other static and dynamic transitions, we consider the 
spin model presented in the next section. 



2. THE SPIN MODEL 

We consider a lattice modelH with, on each lattice site, a Potts variable!^ <Ji, with and 
integer number s of states (a, = l,...,s), coupled to an Ising spin Si with two states 
(Si = ±1). The model can be considered as a schematic representation of structural glasses, 
such as dense molecular glasses, plastic crystal, or orto-therphenyl at low temperature, with 
orientational degrees of freedom frustrated by geometrical hindrance between non-spherical 
molecules. The orientational degrees of freedom are represented by the Potts variables and 
the frustration is modeled by means of ferro/antiferromagnetic interactions for the Ising 
spins. The model is defined by the Hamiltonian 

H s {ri,e i , j } = -2 S jJ2Se ijn ,r j (8) 
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Fig. 1 a) A Potts fully frustrated (PFF) model on a square lattice: on each vertex there is variable 
Ti = SiO~i = ±1, ±2 . . . ± s with s — 4. Here we represent the spin Si by an open or a full dot (respectively 
positive and negative, for example) and the orientation Oi by an arrow pointing in 4 different directions. 
Ferromagnetic (antiferromagnetic) interactions are represented by full (dotted) lines, b) The analytic value 
of the transition temperature (in 2D) T c = T p (q,0), with q = 2s, of the ferromagnetic Potts model (upper 
line), corresponding to X = 0, is given together with the 2D Monte Carlo data for the PSG model (circles), 
corresponding to X = 0.5, with transition at T p PSG and for the PFF model (squares), corresponding to X = 1, 



with transition at Tj, . The data are fitted with Eq.(|lC|) with T* 



The parameters a(X) are shown in the inset. Errors are smaller than symbols size. 



T p ( 9) 0.5) andT;" = T p (q, 1) 



where n = SiOi = ±1, ±2, . . . , ±s has q = 2s states, the sum is extended over all the nearest 
neighbor (n.n.) sites, J is the strength of interaction, eij = ±1 is a quenched variable that 
represents the sign of the ferro/antiferromagnetic interaction, 5 n ^ m = 0, 1 is a Kronecker 
delta. 

In the original formulationu, the Hamiltonian shows a clearer separations between the 
two (cjj and Si) coupled variables: 



tls\Si, &i, €%JJ 



-sJ}^ 5 ai>ff . (ejjSjSj + 1) . 
(hi) 



(9) 



The Eq.(|9|) shows that the Ising spins interactions are diluted by a ferromagnetic Potts 
model, i.e. two n.n. Ising spins can interact only when the corresponding n.n. Potts 
variables are in the same state (or orientation). 

The model has been proposed in two versions. In the first version, the interaction signs 
£j j are randomly assigned, giving rise to frustration and disorder. This model can be 
considered the generalization of an Ising spin glass (SG) - with two states - to a model 
with 2s states - the Potts SG (PSG)i. 

In the other version, the interactions signs are assigned in a deterministic way and there 
is an odd number of €ij = — 1 (antiferromagnetic interaction) on each lattice cell. In this 
way the model has frustration and no disorder and is the generalization of the Ising fully 
frustrated (IFF) modela to a 2s-states Potts fully frustrated (PFF) modelu (Fig.ffl.a). In the 
following we will review some results on the PSG and the PFF model. 
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Fig. 2 The slow relaxation in the PFF model with q = 2s = 4 in 2D. a) The fitting parameter j3 (upper 
panel) and r (lower panel) for the second step of the correlation function of the Potts order parameter 
M (extrapolated to the thermodynamic limit). The circles are the parameters for the Ogielski stretched 
exponential form in Eq.([7J), the triangles for KWW stretched exponential form in Eq.{J6|). The arrow shows 
the estimate of T p (4, 1) in Eq.(|lO|). The lower panel includes also data for the integral correlation time 
T iri t — limt maa .^oo[2 +/Ct=cT /m(£)]. Where not shown, the errors are smaller that the symbols size, b) The 
log- log plot of the autocorrelation time To, for the energy density E (squares) and the Potts order parameter 
M (circles) at the finite-size transition temperature T P (L) as function of the size L = 20, 24, 30, 40, 50. The 
autocorrelation time to is defined as the time (in unit of Monte Carlo steps) at which Jm{jo,T p {L)) = 0.4 
and /_b(to,T p (L)) = 0.3. Here T P (L) — > T p (4, 1) for L — > oo. The positive slopes zm and ze show that these 
times diverge in the thermodynamic limit L — » oo. 

2.1 The Thermodynamic Transitions 

Both the PFF and the PSG model have two thermodynamic transitions, one associated 
to the variables ctj and one to the spins Si. At zero temperature (in 2D) there is a transition 
for the Ising spins S{. The transition is in the IFF universality class for the PFF model, 
and in the Ising SG class for the PSG model. At finite temperature T p it has been shown, 
by means of Monte Carlo (MC) simulations^ and analytic approaches^!, that the Potts 
variables undergo a transition in the universality class of a s-state Potts ferromagnetic 
model (Fig]l|.b). 

It is possible to extend numericallyu the analytic expression of the ferromagnetic Potts 

transition temperaturecJ to the PSG and the PFF case, by introducing the fraction X of 

elementary frustrated cell in the lattice, with the ferromagnetic Potts model corresponding 

to X = 0, the PSG model to X = 0.5 and the PFF model to X = 1. The generalized 

relation is 

k B T p {q,X) 1 



a(X)qJ ln 



\ + ^{X)q 



(10) 



where q = 2s and the parameter a(X) is reported in Fig|l|.b. Note that in 3D the transition 
of the Ising spins is expected at finite T <T p . 



2.2 The Glassy Behavior and the connection with the thermodynamics 

It has been shownHBEl that the PSG and the PFF models have a glassy behavior. For 
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example in the PFF with q = 2s = 2, 4 in 2D and 3D the autocorrelation function 

{A(t,T)A(0,T)) - (A(T)) 2 

IM ' ' (A(0,T) 2 ) - (A(T)) 2 ' { ' 

for a global quantity A (like the total energy E or the Potts order parameters M) shows the 
'two-step' behavior of glasses, with a non-exponential second step (upper panel in Fig.^.a). 

These global correlation functions can be considered as measures of the structural relax- 
ation process. Their relaxation times (tq and Tq 1 in Fig||b) diverge, in the thermodynamic 
limit, at the ordering Potts transition temperature T p (q,l) in Eq.(pI|)El. For T > T p this 
diverging relaxation time resembles the a-relaxation process occurring for T > T g in fragile 
liquids. * 

In the PSG case the onset T* of non-exponentiality coincides with the Griffiths temper- 
ature correspondingI3B to the transition temperature T c (q) = T p (q,0) of the ferromagnetic 
g-states Potts model J This finding generalizes the results for the Ising SGEj'EJ. In the Ising 
SG, indeed, the non-exponential behavior for T < T c was relatecO to the presence of unfrus- 
trated regions, due to the randomness, with a size probability distribution that decreases 
exponentially. At T c the ferromagnetic-like correlation length is equal to the characteristic 
size of each region, giving rise to a non-Gaussian distribution of relaxation times, that is 
responsible for the non-exponentiality of the global relaxation time. 

Due to the lack of randomness, in the PFF model this explanation is not valid (the 
Griffiths temperature is not defined in this case). Indeed, for the PFF model, the onset 
of non-exponential relaxation corresponds to the Potts transition temperature, as has been 
shown numerically&Ei in 2D for q = 2, 4 (upper panel in Fig.|2[a) andH in 3D for q = 2. In 
these cases the previous argument can be extended considering that the ferromagnetic-like 
correlation length is now associated to the Potts variables, ordering at T p . 

Summarizing, the results suggest that if the system is disordered (X ^ 1), then T* 
corresponds to the Potts transition temperature for X = (i.e. the Griffiths temperature 
T c ); if the system is fully frustrated (X = 1), T* corresponds to the Potts transition 
temperature for X = 1 (there is no Griffiths temperature in this case). Note that the 
transition is vanishing for X ^ 1 (because the Griffiths temperature marks a transition 
that disappears for vanishing external field) and for the case X = 1 with q = 2 (in this case 
the transition is defined only as a percolation transition, because there are no orientational 
states). The transition is actually present for the cases with X = 1 and q > 2. 

Hence the onset of the non-exponentiality is marked by a temperature related, at least in 
this model, to a thermodynamic transition. The transition, in this case, is due to the Potts 
(orientational) variables, but, in a more general case, could be associated to any global 

* For T < T p , these correlation times decrease with T (e.g. for M in lower panel in Figfla), because they are 
proportional to £ z , where the correlation length £ is finite for T ^ T p (z is the quantity-dependent dynamical 
exponent, approximated by the exponents in Fig.[2|.b). In real glasses, instead, the relaxation time is always 
increasing for decreasing T. A relaxation time with a more appropriate behavior could be the one associated 
to the overlap between Potts configurations visited at different times. We will propose an 'extension' of the 
model, in the last section, that should give a monotonic slow relaxation time. 

' In disordered systems it is possible to show that for finite external field a free energy (Griffiths) singularity 
ariseaiJ. In the limit of external field going to zero, the temperature at which this singularity occurs goes to 
the transition temperature T c of the corresponding system with no disorder, and the singularity vanishes. 
The case with no disorder, corresponding to the q-states PSG model, is the X — case (i.e. the ferromagnetic 
g-states Potts model) and is T c = T p (q, 0). 
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Fig. 3 The /3-relaxation in the PSG model in 2D for q = 2s — 4. a) The normalized correlation function 
f x (t) of the nonlinear susceptibility: temperatures are in units of J jks ; times in units of MC steps. The 
data are fitted with the forms in Eq.s(H,H,R). The Ogielski form in Eq.(M) is the only able to fit the data 
on 4 decades. Where not shown the error bars are smaller then the symbols size, b) Tha fitting parameters 
{3 and r for the stretched exponential forms in Eq.s(y,[7|). For the values of x see RefJj. The onset T* of 
the stretched exponential corresponds to the Griffiths temperature T c (marked with C) above the Potts 
transition temperature T p (marked with P). The lower panel includes also data for the integral correlation 
time Tint = hm tma!C ^ 00 [i + X^=cT /%(*)]• The lines are only guides for the eyes. 

(structural) ordering. This kind of ordering could be, in principle, due to some internal 
degrees of freedom non easily detected in real experiments. 



3. THE /3-RELAXATION 

Both PSG and PFF models in 2D showo a non-exponential behavior also for the normal- 
ized correlation functions / x (t) of the time-dependent nonlinear susceptibility 



XSG{t) = -jrji 



N 



^2Si(t + to)Si(to) 



%=x 



(12) 



(with N total number of spins, to equilibration time, Xsg(O) = N) where the angular 
brackets stand for the thermal average and the bar stands for the average over the disorder 
in the PSG case, and is absent in the PFF case. Note that f x (t) depends explicitly only on 
the Ising spins. 

In the PSG case, the onset T* of non-exponentiality for fy, corresponds to the Griffiths 
temperature T c (Fig.g). In the PFF case it has been shownQ that T* corresponds to the 
Potts transition temperature. These results seem to have the same interpretation as those 
for the correlation functions of the quantities depending explicitly on the Potts variables 
(previous section). 

The difference in this case is that the Ising spins have no thermodynamic transition at 
T p . Therefore their correlation time r, associated to f x (t), does not diverge at T p (lower 
panel in Fig||.b) and, hence, describes a fast process, respect to the slow dynamics of the 
Potts variables. This fast r can be considered as a measure of the time needed by the Ising 
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Fig. 4 The /3-relaxation in the PFF model in 2D with q = 4. a) The logarithm of the correlation time r, 
associated to f x , for two finite systems with linear sizes L — 40, 120, plotted against the inverse temperature 
(T is in units of J/fcs). The finite size effect is evident only near the Potts transition temperature T p 
(marked by a vertical line) , where the Potts correlation length reach the system size. An Arrhenius behavior 
in this plot is represented by a straight line, b) The data for T > T p plotted against the logarithm of T — T p 
to verify the power law predicted by MCT. The power law is asymptotically satisfied and the slope is an 
approximate estimate of the power exponent. In both panels the lines are only guides for the eyes. 



spins to minimize the energy locally, on the diluted lattice given by the clusters of ordered 
Potts variables. This fast (local) relaxation corresponds to the /^-process. 

The intriguing result in this case is that this correlation time r has a non-Arrhenius 
behavior for T > T p and an Arrhenius behavior for T < T p (Fig|3|.b and Fig^). This 
result resembles the not well understood a — (3 bifurcation, presented in the introduction, 
that is seen at a characteristic temperature T a _p in some fragile liquidsu and intermediate 
liquidsty. 

Moreover the crossover, from non-Arrhenius to Arrhenius, occurs at T p in both models 
(with or without disorder). Therefore T a _p is separated by the onset of non-exponential 
relaxation function T* - marked by T c if there is disorder, and by T p if there is no disorder 
(previous section). 

In the experiments!! the temperature T a _p seems to coincide with the Tmct of the Mode 
Coupling Theory, at which the a-relaxation time diverges as a power law, Eq.(|3|). Here 
the bifurcation coincides numerically with T p , at which the global relaxation time diverges. 
Indeed, the correlation time r can be fitted with the VTF law, Eq.(J2J), with To = T p , or can 
be (asymptotically) well described by a power law diverging at T p (e.g. Fig.||.b). 

These results suggest that the bifurcation here is related to the thermodynamic transition 
of the Potts variables. The fast /3-process is driven by the slow dynamics of the global process 
for T approaching T p from above. The non-Arrhenius behavior of the global process induces 
a non-Arrhenius behavior in the coupled fast variable and the dominant dynamic process 
is the one related to the slow variable. As has been shown for supercooled liquidsEJ, is 
reasonable that in this regime the dominant dynamic process is the relaxation to the lower 
accessible energy, at the given T, following the instability directions, in the energy landscape, 
of the visited states, i.e. the activated processes play no role. 

For the supercooled liquid caseO it has been shown that Tmct corresponds to the temper- 
ature below which the average number of instability directions for the visited configuration 
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goes to zero. From preliminary results, the same situation appears to be valid here at T p . 

Finally, below T p the system is exploring one of the basins of attraction of the Potts 
variables a% and the dynamics of the fast Ising spins is no longer coupled to that of the Ui. 
The dynamics of the fast variables become activated (Arrhenius) consistently with the lack 
of instability directions. Indeed, in this case the energy cannot be lowered by following an 
instability direction, but only by doing an activated process. This is analogous of what has 
been found for the global relaxation in supercooled liquids", but in our case the focus is 
on the local /^-processes. 

Therefore the results here suggest that the a — (3 bifurcation can be related to the ther- 
modynamics of the systems, and that the Potts transition temperature T p of these systems 
plays the role of the Tmct of supercooled liquids. 



4. CONCLUSION AND PERSPECTIVES 

The glassy spin models considered here, the PSG and the PFF model, are characterized 
by the presence of two coupled variables, whose dynamics decouple when the relaxation 
time of the slow variable diverges. This happens at the temperature T p where the slow 
variable has a thermodynamic phase transition. 

The picture that can be derived is the following. Due to the symmetry breaking of the 
slow variable for T < T p , the system is attracted in one of the basin of the energy landscape 
of the slow variable. Inside this basin, the dynamics of the system is mainly due to the fast 
variable and the dominant dynamical mechanism is the activated process, as consequence 
of the lack of instability directions for the visited states. 

Macroscopically, this is shown by the a-/3 bifurcation that can be seen in some fragile 
and intermediate glassforming liquids&EI, i.e. by the crossover of the /3-relaxation time from 
a non- Arrhenius behavior to an Arrhenius law at a temperature T a _^. 

The experiments suggest that T a _p ~ Tmct, consistently with the analysis reported 
here. Indeed, in our models the bifurcation occurs at the temperature which plays the role 
of Tmct, i-e. T p , where the slow process has a diverging relaxation time. 

Therefore the a-(3 bifurcation can been related to the thermodynamics, as well as the 
onset T* of non-exponential relaxation. What we can learn from the 'toy model' studied 
here is that, in cases in which there is disorder, T* and T a _p do not coincide. The first, 
indeed, corresponds to the Griffiths temperature T c - associated to a real transition only for 
a non-zero external field coupled to the slow variable - and the second to the temperature 
T p , where the slow relaxation time diverges. The two temperatures coincides only in the 
particular case of a frustrated system with no disorder. 

An open questions is the effect of the T-dependence of the global correlation length 
£. An interesting case is the one in which £ does not decrease below the thermodynamic 
transition temperature. A model with such characteristic is the XY model, undergoing the 
Kosterlitz-Thouless-BerezinskiiEO transition. The resulting Hamiltonian will be 

H{Si,4>i,eij} = -J^cos(0j - 4>j - Aij)(eijSiSj + 1) (13) 

where Aij are constants depending on the gauge (or the external field) and fa £ [0, 2n) are 
a more realistic representation of the continuous orientational variables. 
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